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ABSTRACT 

The magnetic field in the solar convective zone has a random small-scale com¬ 
ponent with the r.m.s. value substantially exceeding the strength of a regular 
large-scale field. For two Majorana neutrino flavors x two helicities in the pres¬ 
ence of a neutrino transition magnetic moment and nonzero neutrino mixing 
we analyze the displacement of the allowed (Am^ -sin^ 9) - parameter region 
reconciled for the SuperKamiokande(SK) and radiochemical (GALLEX, SAGE, 
Homestake) experiments in dependence on the r.m.s. magnetic field value b. We 
find an effective production of electron antineutrinos in the Sun even for small 
neutrino mixing through the cascade conversions in the random magnetic field 
that would be a signature of the Majorana nature of neutrino if hen will be 
registered. Basing on the present SK bound on electron antineutrinos we also 
estimated the corresponding excluded area in the plane Amf, sin^ 29. 

1. Introduction and stating the problem 

Recent results of the SuperKamiokande (SK) experiment0 have already confirmed 
the solar neutrino dehcit at the level less Rsk/ Rssm ^ 0.4. 

Moreover, day/night (D/N) and season neutrino flux variations were analized in 
this experiment and in 4 years one expects to reach enough statistics in order 
to conhrm or to refuse these periods predicted in some theoretical models. D/N 
variations would be a signature of the MSW solution to the Solar Neutrino Problem 
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(SNP). The signature of the vacuum oscillations is a seasonal variation of neutrino 
flux, in addition to geometrical seasonal variation l/r^(f). And signature of the 
resonant spin-flavor precession (RSFP)-solution is the 11 year periodicity of neutrino 
flux. 

In all four solar neutrino experiments one measures the integral spectrum, i.e. the 
number of neutrino events per day for the SK experiment, 

N" = EE# / UE)aa{E) < Paa{E) > dE , (1) 

a i JEthiSK) 


or the number of neutrino events in SNU (1 SNU = 10 captures per atom/per sec) 
in the case of GALLEX (SAGE) and Homestake experiments. 
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Ai(E)crf“’ < PUE) > dE , 


( 2 ) 


where is the integral flux of neutrinos of kind ”i” {i = pp , Be ,pep , B) assumed 
to be constant and uniform at a given distance r = t from the center of the Sun, 
Xi{E) is the normalized differential flux, Xi(E)dE = 1 ; Ca are corresponding 

cross sections; the thresholds Eth{Ga,Cl) for GALLEX (SAGE) and Homestake are 
0.233 MeV and 0.8 MeV correspondingly while for SK at present Tth = 6.5 MeV i. 
The probabilities Paaif) = where the subscript a equals to a = e for z/eu, 

a = e for PeRi a = /i for and a = /I for correspondingly, satisfy the unitarity 
condition 

Pee{t) + Pee{t) + p Pfip,{t) = 1 , (3) 

and “< ... >” corresponds to averaging over the solar interior if necessary. 

All three elementary particle physics solutions (vacuum, MSW and RSFP) suc¬ 
cessfully describe the results of SK, GALLEX, SAGE and Homestake solar-neutrino 
experiments, because the suppression factors of neutrino fluxes are energy dependent 

with high precision 

at all radial distances of interest (see recent review by Berezinsky i). 

There exists, however, a problem with large regular solar magnetic helds, essential 
for the RSFP scenario. It is commonly accepted that magnetic helds measured at the 
surface of the Sun are weaker than within interior of the convective zone where this 
held is supposed to be generated. The mean held value over the solar disc is about 
of order 1 G and in the solar spots magnetic held strength reaches 1 kG. 

Because sunspots are considered to be produced from magnetic tubes transported 
to the solar surface due to the boyancy, this hgure can be considered as a reasonable 
order-of-magnitude observational estimate for the mean magnetic held stren^h in 
the region of magnetic held generation. In the solar magnetohydrodynamics □ one 
can explain such helds in a self-consistent way if these helds are generated by dynamo 
mechanism at the bottom of the convective zone (or, more specihc, in the overshoot 
layer). But its value seems to be too low for ehective neutrino conversions. 


and helioseismic data conhrm the Standard Solar Model (SSM) 
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The mean magnetic field is however followed by a small scale, random magnetic 
field. This random magnetic field is not directly traced by sunspots or other tracers 
of solar activity. This field propagates through convective zone and photosphere 
drastically decreasing in the strength value with an increase of the scale. According 
to the available understanding of solar dynamo, the strength of the random magnetic 
field inside the convective zone is larger than the mean field strength. A direct 
observational estimation of the ratio between this strengthes is not available, however 
the ratio of order 50 - 100 does not seem impossible. At least, the ratio between the 
mean magnetic field strength and the fluctuation at the solar surface is estimated as 
50 (see e.g. H). 

This is the main reason why we consider here an analogous to the RSFP scenario, 
an aperiodic spin-flavour conversion (ASFC), based on the presence of random mag¬ 
netic fields in the solar convective zone. It turns out that the ASFC is an additional 
probable way to describe the solar neutrino deficit in different energy regions, espe¬ 
cially if current and future experiments will detect electron antineutrinos from the 
Sun leading to conclusion that neutrinos are Majorana particles. The termin “ape¬ 
riodic” simply reflects the exponential behaviour of conversion probabilities in noisy 
media (cf. 1,0). 

As well as for the RSFP mechanism all arguments for and against the ASFC 
mechanism with random magnetic fields remain the same ones that have been recently 
summarized and commented by Akhmedov (see i and references therein). 

But contrary to the case of regular magnetic fields we find out that one of the 
signatures of random magnetic field in the Sun is the prediction of a wider allowed 
region for neutrino mixing angle in the presence of from the Sun including the case 
of small mixing angles. 

2. Magnetic Field Model and Numerical Simulation 

The random magnetic field is considered to be maximal somewhere at the bottom 
of convective zone and decaying to the solar surface. To take into account a possibility, 
that the solar dynamo action is possible also just below the bottom of the convective 
zone (see 0), we accept, rather arbitrary, that it is distributed at the radial range 
0.7^0 — I.ORq, i.e. it has the same thickness as the convective zone, while its 
correlation length is Lq = 10^ km that is close to the mesogranule size. 

We suppose also that all the volume of the convective zone is covered by a net of 
rectangular domains where the random magnetic field strength vector is constant. The 
magnetic field strength changes smoothly at the boundaries between the neighbour 
domains obeying the Maxwell equations. Since one can not expect the strong influence 
of small details in the random magnetic field within and near thin boundary layers 
between the domains this oversimplified model looks applicable. 

In agreement with the SSM i the neutrino source is supposed to be located inside 
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Fig 1. Geometry of neutrino trajectories in random magnetic fields. Finite 
radius of r'eL-neutrino source, the solar core radius O.li?©, is shown. 

the solar core with the radius of the order of Ry = O.li?© = 7 x 10"^ km. For the sake 
of simplicity we do not consider here the spatial distribution of neutrino emissivity 
from a unit of a solid angle of the core image and the differences in for different 
neutrino kinds ”i”. Different parallel trajectories directed to the Earth cross different 
magnetic domains because the domain size Lq (in the plane which is perpendicular 
to neutrino trajectories) is much less than the transversal size = Ry of the full set 
of parallel rays, Lq -C Ry, see Fig. 1. The whole number of trajectories (rays) with 
statisticaly independent magnetic helds is about R^/Lq ~ 50. 

At the stage of numerical simulation of the random magnetic field we generate the 
set of random numbers with a given r.m.s. value and for each realization of random 
magnetic fields along each of 50 rays solve the Cauchi problem for the (4 x 4) master 
equation , 



^ keL ^ 
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where C 2 = cos26', S 2 = sin 20, 5 = Am‘^/4E are the neutrino mixing parameters; /i = 
p .12 is the neutrino active-active transition magnetic moment; H±{t) = Bo±it) + b±{t), 
H± = Hx ± iHy, are the magnetic held components consisting of regular {Bo±) and 
random {b±) parts which are perpendicular to the neutrino trajectory in the Sun; 
Ve{t) = GFV2{p{t)/mp){Ye - Yn/2) and V^{t) = GFV2{p{t)/mp){-Yn/2) are the 
neutrino vector potentials for Vei and in the Sun given by the abundances of the 
electron (17 = '^pNeit)/p{t)) and neutron (l^ = mpNn{t) /p(t)) components and by 
the SSM density prohle p{t) = 250 g/cm^ exp(—10.54t) 0. Then we calculate the 
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dynamical probabilities Paa{t), and average them in transversal plane, i.e. obtain the 
mean arithmetic probabilities as fnnctions of mixing parameters sin^ 26 , 6. We argue 
that being additive functions of the area of the convective zone layer (or, the same, 
of the number of rays) these probabilities with increasing area become self-averaged. 


3. Asymptotic solution 


To demonstrate the properties of our model we briefly discuss how random mag¬ 
netic held inhuences the small-mixing MSW solution to SNP. It turns out that for the 
SSM exponential density prohle, typical borone neutrino energies ~ 7 lAMeV, 
and AnP ~ 10“^ eV'^, the MSW resonance (i.e. the point where = 2c2S) 

occurs well below the bottom of the convective zone. Thus we can divide the neutrino 
propagation problem and consider two successive stages. First, after generation in 
the middle of the Sun, neutrinos propagate in the absence of any magnetic helds, 
undergo the non-adiabatic (non-complete) MSW conversion and acquire 

certain nonzero values Pee end P^^, which can be treated as initial conditions at the 
bottom of the convective zone. For small neutrino mixing S 2 ^ 0 the (4 x 4)-master 
equation Eq.( then splits into two pairs of independent equations describing cor¬ 
respondingly the spin-havor dynamics UeL and —>• he_R in noisy magnetic 

helds. In addition, once the MSW resonance point is far away from the convective 
zone, one can also omit 14 and I 4 in comparison with C 2 ^. For ^iiR conversion 

this results into a two-component Schrodinger equation 



-6 

/i6+(f) 
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5 
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Pee(O), 


= Pfifi{0) 

= 0 . As normalized 


with initial conditions |z/eLp( 0 ) = 
probabilities Peeif) and Pfifiit) (satisfying the conservation law Peeif) + Pfifiit) = 
Pee(O)) are the only observables, it is convenient to recast Eq.(|^) into an equivalent 
integral form 


s{t) = S{0)-4fi‘^JdpJdt2S{t2) 

{[bx{ti)K{t2) + by{ti)by{t 2 )] cos25(4 - ^ 2 ) 

-[bx{ti)byit 2 ) - bxit 2 )by{ti)] sin25(fi - 4)}, 

where S(t) = 2Pee(t) — Pee(O) = Peeit) “ Pfifiit) is the third component of the polar¬ 
ization vector S. 

Let us trace the line of further derivation in a spirit of our numerical method. 
Dividing the interval of integration into a set of equal intervals of correlation length 
Lq, and assuming that possible correlations between S{t 2 ) and 4, bj under the integral 
are small, S{t) itself varies very slowly within one correlation cell, and making use 
of statistical properties of random helds, {i) diherent transversal components within 
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one cell are independent random variables, and {ii) magnetic fields in different cells 
do not correlate, we can average Eq.(^ thus obtaining a finite difference analogue. 
Returning back to continuous version we get 

S(t) = S(0)exp j < Pit') > dA , (7) 

where we retained possible slow space dependence of the r.m.s. magnetic field value. 

For S <C Rq ^ constant r.m.s. < 6^ > we obtain the simple ^-correlation re¬ 
sult 1,0, otherwise there remains an additional stabilizing factor sin^a:/x^, demon¬ 
strating to what extent the vacuum/medium neutrino oscillations within one corre¬ 
lation cell can suppress the spin-flavor dynamics due to the magnetic field only. 

Another important issue is the problem of temporal dependence of higher statis¬ 
tical moments of P^a- As Paa enter the Eq.(|^ for the number of events one should be 
certain that the averaging procedure does not input large statistical errors, otherwise 
there will be no room for the solar neutrino puzzle itself. For 6 -C basing on Eq. 
( we evaluate the dispersion cXp. Here we present only the final result, the details 
of derivation will be published elsewhere: 

<Tp = Pl(t) >-< P„(t) >= = ^ (l - , (8) 


< > - ~2^ (l + e , (9) 

where T = (4/3)/i^ < > Lq. We see that relative mean square deviation of Pee 

from its mean value tends with t —>■ cx) to its maximum asymptotic value 


(Tpjt) 

< Pee{t) > 


0.707, 


f —>• CXO, 


( 10 ) 


irrespectively of the initial value Pee(O). 

This estimate is true, evidently, only for one neutrino ray. Averaging over N 
independent rays lowers the value (p!o|) in a/]V times. That is for our case of ~ 50 
rays we get that maximum relative error should not exceed approximately 10%, thus 
justifying the validity of our approach. For smaller magnetic fields the situation is 
always better. 

To conclude this section, it is neccesary to repeat that the above estimate Eq. 
([T0|) indicates possible danger when treating numerically the neutrino propagation 
in noisy media. Indeed, usually adopted one-dimensional (i.e. along one ray only) 
approximation for the (4 x 4) master equation Eq. or (2 x 2) Eq. (|]) can suffer 
from large dispersion errors and one should make certain precautions when averaging 
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these equations over the random noise before numerical simulations. Otherwise, the 
resulting error might be even unpredictable. 


4. Discussion of results 


In order to find the regions in the Amf, sin^ 29 plane excluded from nonobservation 
of UeR in SK it is necessary to compute and plot the isolines of the ratio . 

Indeed, as antineutrinos were not detected in the SK experiment one should claim 
that antineutrino fiux is smaller than the SK background > 8-3 MeV) < 

6 X 10^ i. Deviding this inequality by the SSM borone neutrino fiux with 

energies E > Eth = 8.3 MeV, > 8.3 MeV) = 1.7 x 10® we 

find the bound on the averaged over cross-section and spectrum, cascade transition 
probability VeL i>eR , 


^pXE > 8.3 MeV) _ fll\B{E)P,,{E)ap{E)dE 
^f^{E > 8.3) fll XB{E)ap{E)dE 




(a) 


(b) 


Fig 2. The averaged (over 50 rays) probability Pgg = i^eR^eR-) ™ dependence 
on 5 = Am?/AE and sin^ 20-parameters for br.m.s. = 50 kG (a) and hr.m.s. = 
100 kG (b). The dimensionless 6 = 8.81 x Am?{eV‘^)/E{MeV) is plotted 
along y-axis. 
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or < 3.5% . Here au{E) = 9.2 x - 1.3 MeH)/10 MeVf is 

the cross-section of the capture reaction Vf.p —> ne^. Notice that the above bound 
is not valid for low energy region below the threshold of the antineutrino capture by 
protons, E < 1.8 MeV. 

For low magnetic fields, B ,b < 20 kG, both for regular and random ones we 
do not find any violation of the bound Eq. (Q. However, for strong magnetic 
helds such forbidden parameter regions appear in different areas over 6 and sin^ 29 
for different kind of magnetic helds (regular and random). In general, the inhuence 
of random magnetic helds is more pronounced as compare to the regular ones of the 
same strength. Here we illustrate only the random magnetic held case, see Fig. 2,3 
for brms = 50 and 100 kG. It follows that the more intensive the r.m.s. held ( 6 ^) in 
the convective zone the more ehective spin-havour conversions lead to the production 
of the right-handed UeR, t'^/j-antineutrinos. In the case of small mixing this is also 
seen from Eq. (^. 
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Fig 3. The reconcilement of all experiments in the plane Am^(eF^) and sin^ 29. 
The isolines with dashed bands represent the DATA/SSM ratio equal corre¬ 
spondingly to 0.368 ± 0.026 (SK), 0.509 ± 0.059, 0.504 ± 0.085 (GALLEX 
SAGE), and 0.274 ± 0.027 (Homestake). 

(a) The random field brms = 50 kG; (b) the random field brms = 100 kG. 






















In Fig. 3 we present the allowed parameter regions to reconcile four solar neutrino 
experiments along with the SK bound on Vf.R given by Eq. O- One can see that 
without the latter bound there exist two commonly adopted parameter regions, the 
small mixing and the large mixing ones, currently allowed or excluded in dependence 
of the values of the r.m.s. magnetic field parameters. Here we demonstrate the allowed 
case only, when excluded by the SK bound region is well below the two parameter 
regions of interest. 

Our results, however, show that there exists strong dependence not only on the 
r.m.s. magnetic field strength but on its correlation length also. In particular, when 
for the same brms = 100 kG we substitute Lq = 10^ km (a granule size) instead of 
Lq = 10^ km (a mesogranule size treated in numerical simulation) the small mixing 
MSW region is fully excluded via the SK bound Eq. (|TI[). This is easily explained 
due to the maximum value of ASFC seen from Eq. (sinx/x —> 1). Vice versa, if 

the small-mixing MSW solution is valid and a large neutrino magnetic moment exists 
we can extract from neutrino data important constraints on the structure of random 
magnetic helds deep under surface of the Sun. 


5. Conclusions 

We develop a model of neutrino spin-flavour conversions in the random magnetic 
fields of the solar convective zone supposing in consistence with modern MHD models 
of solar magnetic fields that random fields are naturally much higher than large-scale 
magnetic fields created and supported continuously from the small-scale random 
ones 0. 

It follows that if neutrinos have a large transition magnetic moment their dynam¬ 
ics in the Sun is governed by random magnetic helds that , hrst, lead to aperiodic 
and rather non-resonant neutrino spin-havor conversions, and second, inevitably lead 
to production of electron antineutrinos for low energy or large mass difference re¬ 
gion. If antineutrinos UeR would be found with the positive signal in the BOREXINO 
experiment^ or, in other words, a small-mixing MSW solution to SNP fails, this will 
be a strong argument in favour of magnetic held scenario with ASFC in the presence 
of a large neutrino transition moment, /i ~ for the same small mixing angle. 

The search of bounds on /i at the level /i ~ in low energy z/ e-scattering, 

currently planning in laboratory experiments 0, will be crucial for the model consid¬ 
ered here. 

We would like to emphasize the importance of future low-energy neutrino ex¬ 
periments (BOREXINO, HELLAZ) which will be sensitive both to check the MSW 
scenario and the z/gij-production through ASFP. As it was shown in a recent work0 
a diherent slope of energy spectrum prohles for diherent scenarios would be a crucial 
test in favour of the very mechanism providing the solution to SNP. 
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